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Abstract. We present the first results on galaxy metallicity evolution at z>3 from two projects, 
LSD (Lyman-break galaxies Stellar populations and Dynamics) and AMAZE (Assessing the 
Mass Abundance redshift Evolution). These projects use deep near- infrared spectroscopic ob- 
servations of a sample of ~40 LBGs to estimate the gas-phase metallicity from the emission 
lines. We derive the mass-metallicity relation at z>3 and compare it with the same relation at 
lower redshift. Strong evolution from z=0 and z=2 to z=3 is observed, and this finding puts 
strong constrains on the models of galaxy evolution. These preliminary results show that the 
effective oxygen yields does not increase with stellar mass, implying that the simple outflow 
model does not apply at z>3. 
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1. Introduction 

Metallicity is one the most important property of galaxies, and its study is able to shed 
light on the detail properties of galaxy formation. It is an integrated property, not related 
to the present level of star formation the galaxy, but rather to the whole past history of 
the galaxy. In particular, metallicity is sensitive to the fraction of baryonic mass already 
converted into stars, i.e., to the evolutive stage of the galaxy. Also, metallicity is effected 
by the presence of inflows and outflows, i.e., by feedback processes and the interplay 
between the forming galaxy and the intergalactic medium. 

It is well known that local galaxies follow a well-defined mass-metallicity relation, 
where galaxies with larger stellar mass have higher metallicities (jTremonti et al. 2 004, 
|Lee et al., 2006[ ). The origin of the relation is uncertain because several effects can be, 
and probably are, active. It is well known that in the local universe starburst galaxies 
eject a significant fraction of metal-enriched gas into the intergalactic medium because 
of the energetic feedback from exploding SNe, both core-collapse and, possibly, type la 
(jMannucci et al. 2006^ . Outflows are expected to be more important in low-mass galaxies, 
where the gravitational potential is lower and a smaller fraction of gas is retained. As a 
consequence, higher mass galaxies are expected to be more metal rich (see, for example, 
Ed munds 19901 [Garnett 20 02). A second possibility is related to the well know effect of 
"downsizing" (e.g.. lCowie et al. 19961) . i.e., lower-mass galaxies form their stars later and 
on longer time scales. At a given time, lower mass galaxies have formed a smaller fraction 
of their stars, therefore are expected to show lower metallicities. Other possibilities exist, 
for example some properties of star formation, as the initial mass function (IMF), could 



change systematically with galaxy mass (Koppen et al. 2007). 
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All these effects have a deep impact on galaxy formation, and the knowledge of their 
relative contributions is of crucial importance. Different models have been built to re- 
produce the shape of the mass-metallicity relation in the local universe, and different 
assumptions produce divergent predictions at high redshifts (z>2). To explore this issue 
several groups have observed the mass-metallicity relation in the distant universe, around 
z=0.7 dSavaglio et al. 2005[ ) and z=2.2 (|Erb et al. 2006|1 . They have found a clear evolu- 
tion with cosmic time, with metallicity for a given stellar mass decreasing with increasing 
redshift. 

For several reasons, it is very interesting to explore even higher redshifts. The red- 
shift range at z^3-4 is particularly interesting: it is before the peak of the cosmic star 
formation density (see, for example. IMannucci et al. 2007[) , only a small fraction ( 15%, 
IPozzetti et al. 2007[) of the total stellar mass have already been created, the number of 
mergers among the galaxies is much larger than at later times (jConselice et al. 2007[) . 
As a consequence, the prediction of the different models tend to diverge above z=3, and 
it is important to sample this redshift range observationally. The observations are really 
challenging because of the faintness of the targets and the precision required to obtain 
a reliable metallicity. Nevertheless, the new integral-field unit (IFU) instruments on 8-m 
class telescopes are sensitive enough to allow for the project. 



2. LSD and AMAZE 

Metallicity at z~3 can be obtained by measuring the fluxes of the main optical emis- 
sion lines ([Oil], H/3, [OIII], Ha), whose ratios have been calibrated against metallicity 
in the local universe ( |Nagao et al. 200 6, Kewley & Ellison 2008). Of course, this method 
can be applied only to line-emitting galaxies, i. e., to low-extinction, star- forming galax- 
ies, whose line can be seen even at high-redshifts. In contrast, the gas metallicity of 
more quiescent and/or dust extincted galaxies, like EROs (jMannucci et al. 2002[) , DRGs 
(jFranx et al. 2003)) . and SMGs (Chapman et al. 2005), cannot be easily measured at high 
redshifts Q- 

For the AMAZE project we observed 30 galaxies at z^3.3 from various sources (see 
IMaiolino et al. 2008l for details). We only selected galaxies having a good SPITZER/IRAC 
photometry (3.6-8 /im), an important piece of information to derive a reliable stellar 
mass. These galaxies were observed, in seeing-limited mode, with the integral-field unit 
(IFU) spectrometer SINFONI on ESO/VLT, with integration times between 3 and 6 
hours. When computing line ratios, it is important that all the lines are extracted from 
exactly the same aperture, to avoid differential slit losses that could spoil the line ratio. 
In this, the use of an IFU is of great help. We observed the H and K bands simultaneously 
with spectral resolution R^1500, providing a full coverage of all the most important lines. 
This paper is based on about 1/3 of the full data sample, while the remaining fraction is 
still under analysis. 

For LSD, we extracted 10 galaxies from the UV-selected Lyman Break Galaxies (LBG) 
sample by Stei del et al. 20041 The aim of this project is not only to measure metallicities, 
but also to obtain spatially-resolved spectra to measure dynamics and spectral gradients. 
For this reason we used adaptive optics to obtain diffraction-limited images and spectra 
in the near-IR. The target galaxies were chosen to be within 30" of bright foreground 
stars, needed to drive the adaptive-optics system. As the presence of a nearby bright 
star is the only request, this sample, albeit small, is expected to be representative of 



f Stellar metallicities, measured by absorption lines, can also be measured if enough observing 
time is provided, and will be subject of a future work. 
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Figure 1. Composite spectra from AMAZE (top) and LSD (bottom). 

the full population of the LBGs. For LSD we used SINFONI with the same resolution of 
AMAZE and similar integration times. About half of the LSD galaxies have been already 
analyzed. The typical spatial resolutions is about 0.2". 

Fig. [2 shows the composite spectra of AMAZE (top) and LSD (bottom). For both 
AMAZE and LSD, the main optical lines are detected in most of the galaxies. A few 
targets also show the [NeIII]3869 line, and important piece of information to derive a 
robust measurement of metallicity ( |Nagao et al. 2 006). 



3. The Mass- metallicity relation and the effective yields 

Stellar masses are derived by fitting the spectral-energy distributions (SEDs) with 
spectrophotometric models of galaxy evolutions, as detailed in IPozzetti et al. 2 007 and 
IGrazian et al. 20071 These fits also provide estimates of the age and dust extinction of 
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the dominant stellar population. The presence, for most of the objects, of good IRAC 
photometry, allows the determination of reliable stellar masses as the SED is sampled up 
to the rest-frame J band. 

Metallicities are derived by a simultaneous fit of all the available line ratios, as ex- 
plained Maiol ino et al. 2 008. In practice, the derived value is determined by the R23 
indicator or, similarly, by the [OIIIJ5007/H/3 ratio, while the [OIII] 5007/ [Oil] 3727 ratio is 
used to discriminate between the two possible branches of these ratios. The [NeIII]3869/[OII]3727 
line ratio is also very important, when both lines are detected. The uncertainties on metal- 
licity are due to both the spread of the calibration and on the observational error on the 
line ratio, and on average amount to 0.2-0.3 dex. 

Fig. [2] shows the resulting mass-metallicity relation, compared to the same relation as 
measured at lower redshift. A quantitative interpretation of this result will be given when 
the whole data sample will be analyzed. Nevertheless, a strong metallicity evolution can 
be seen, i.e., galaxies at z^3.1 have metallicities ~6 times lower than galaxies of similar 
stellar mass in the local universe. 

The presence, at z>3, of galaxies with relatively high stellar masses (log(M/M©)=9- 
11) and low metallicity put strong constraints on the process dominating galaxy for- 
mation. Several published models fe.g.. lde Rossi et al. 20071 |Kobayashi et al. 2007[ ) can- 
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Figur e 2. Evolution of th e mass- metallicity relatio n from z=0.07 ( |Kewley fc E llison 2 008[ ), to 
z=0.7 dSavaglio et al. 2005| ), z=2.2 l|Erb et al. 2006)1 and z=3.1 (AMAZE+LSD). All data have 
been calibrated to the same scale in order to make all the different results directly comparable. 
Crosses show the AMAZE galaxies, empty dots the LSD galaxies. Solid dots with error bars show 
the average metallicity of the galaxies in each stellar-mass bin, with the associated dispersion. 
The lines show quadratic fits to the data. The line corresponding to z=3.1 shows the fit in 
Maiolino et al. (2008). 
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not account for such a strong evolution, and the physical reason for this can be traced 
to be due to some inappropriate assumption, for example about feedback processes or 
merging history. When taken at face value, some other models (e.g., IBrooks et ai~2 007. 
ITornatore et al. 2007[) provide a better match with the observations, but a meaningful 
comparison can only be obtained by taking into account all the selection effects and 
observational biases. 

Important hints on the physical processes shaping the mass-metallicity relation can be 
understood by considering the effective yields, i.e., the amount of metals produced and 
retained in the galaxy per unit mass of formed stars. If outflows are the main effect in 
shaping the mass-metallicity relation, then the effective yields are expected to increase 
with stellar mass because all the galaxies have converted, on average, the same fraction 
of baryonic mass into stars, but lower mass galaxies have lost a larger fraction of metals 
into the intergalactic medium. This is what is observed in the local universe by Tremonti 
et al. (2004) (but not by |Lee et al., 2006[ ). For LSD and AMAZE, preliminary results 
show that the effective yields tend to decrease, rather than increase, with stellar mass, 
as in Erb et al. (2006). Such a decrease is partly due to selection effects (see Mannucci et 
al., in preparation, for a full explanation), but the observed trend seems to be larger than 
what can be attributed to the effects of biases. The presence of higher yields at lower 
stellar mass imply that outflows decreasing with galaxy mass are not the main driver of 
the mass-metallicity relation at z>3, and different possibilities, related to the efficiency 
of star formation fe.g.. IBrooks et aL~2 007) or different outflow schemes (|Erb et al. 2006[) 
must be considered. 
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